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INTRODUCTION 

Many future NASA missions would utilize significantly 
large and flexible spacecrafts and would require very 
stringent pointing and vibration suppression requirements . 
The active controller that can achieve these objectives will 
have to be designed with very accurate knowledge of the 
dynamic behavior of the spacecraft to ensure performance 
robustness to a variety of disturbances and uncertainties. 

In the past few years, several design approaches were 
proposed for vibration cr-'trol during and after slewing 
maneuvers. NASA Langley Research Center initiated the 
Spacecraft Control Laboratory Experiment (SCOLE) program [1] 
to promote direct comparision and realistic test of various 
control design techniques against a common laboratory 
article. The article was intended to resemble a large space 
antenna attached to the space shuttle orbiter by a long 
flexible mast . 

The primary control objective of SCOLE is to direct the 
RF 1 ine-of-sight (LOS) of the antenna-like conf iguration 
towards a fixed target under the conditions of minimum time 
and limited control authority. 

This problem of directing the LOS of antenna-like 
configuration is studied as being composed of two control 
phases during this research period. In the first phase, the 
LOS of antenna-like structure is to be changed according to 
a prespecified target (slewing) together with minimization 
of vibration amplitudes to ensure stability. The second 
phase is to achieve the total vibration suppression at the 
end of slew maneuver with an augmented control law. ' 

BRIEF SUMMARY OF MAJOR ACCOMPLISHMENTS 

(a) Nonlinear Model Development: 

The focus of the earlier part of the research was in 
upgrading the dynamics of the SCOLE model to reflect all the 
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kinematic nonlinearities in the previous dynamical model. 

The dynamical equations of slewing maneuvers of this large 
flexible spacecraft were developed by formulating Lagranga'n 
equations using an 4 nertial co-ordinate system and a body- 
fixed coordinate system at the point of attachment of the 
flexible beam to the shuttle. The generic model used for 
this analysis consisted of distributed parameter beam with 
two end masses. The three dimensional linear vibration 
analysis of this free-frse beam model with end masses [2] 
was incorporated together with rigid-body slewing maneuver 
dynamics [3] to yield the final set of highly nonlinear and 
coupled equations. 

(b) Slew Maneuvers 

The slew maneuvers were analyzed using this enhanced 
dynamical model in terms of both pure rigid-body slew 
maneuvers and rigid-body slew maneuvers together with the 
suppression of the first two flexible modes. The dynamics of 
the motion during this phase wa$ derived in terms of four 
Euler parameters, and using the method of nonlinear 
decoupling, this nonlinear slew maneuver control problem was 
reduced to calculating a pair of constants for each of the 
first three Euler parameters and the flexible modes in 
implementing output feedback. 

(c) Vibration Suppression at the End of Slew Maneuver 

Since the oscillations of the beam at the end of the 
slew maneuver were considered to be small, the problem of 
vibration suppression was considered to be a linear control 
problem. Thus it was formulated as an infinite-time 
regulator problem with control spillover and observation 
spillover terms incorporated into the performance index 
term. The resultant optimal state feedback control law 
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achieved total vibration suppression. However, the coupling 
among the inodes was found to be a significant problem. 

A paper describing some of the details is included in 
the appendix. This paper is to be presented at the AIAA 
Guidence, Navigation and Control Conference in Williamsburg, 
Va. August 18-20, 1986 and is to be published in the volume. 
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Abstract: 

In this paper, the dynamics and control of slewing 
maneuvers of a large flexible spacecraft namely, NASA- 
Spacecraft Control Laboratory Experiment (SCOLE) test 
article are studied. The dynamical equations obtained 
for slewing maneuvers are highly nonlinear and coupled. 
The maneuver is expressed in terms of four Euler 
parameters and is specified as the angular displacement 
about an arbitrary axis. The slew maneuver control 
problem is developed in terms of rigid-body slewing and 
suppression of two elastic modes is analyzed using 
the method of nonlinear decoupling. 

Nomenclature: 

a(z) - Position vector of masselement on the beam 

— from the point of attachment 

q - Position vector from the point of attachment 

“ to the mass center 

d(z,t) - Displacement vector of mass element in the 
body- fixed frame 

Fg (t } - Force applied at the orbiter mass center 
F (t) - Force applied at the reflector mass center 
Go(t) - Moment applied about the orbiter mass center 

11 - Mass moment of inertia matrix of the shuttle 

1 2 - Mass moment of inertia matrix of the 

reflector 

L - The length of the beam 

m - Total mass of the flexible beam 

mi - Mass of the orbiter 

m2 - Mass of the reflector 

R - Position vector of the mass center of the 

orbiter in the inertial frame 
t - Position vector from the orbiter 

mass center to the point of attachment 
u Jz,t)- The beam deflection in x direction referred 
to the body-fixed frame 

u fz,t)- The beam deflection in y direction referred 
^ to the body-fixed frame 

u Q (z,t)- The torsional deflection about z axis in the 
body-fixed frame 

V - Velocity vector of the mass center of the 

orbiter in the body-fixed frame 
Vq - Velocity vector of the point of attachment 

— in the body-fixed frame 

p - Mass per unit length 

; - Vector representing the axis of rotation 




^ - The angular* velocity of the orbiter in the 

body-fixed frame 
ip - Angle of rotation 
- Introduction; 

Many future NASA missions would utilize 
significantly large and flexible spacecrafts and would 
require very stringent pointing and vibration 

suppression requirements. The active controller that 
can achieve these objectives will have to be designed 
with very accurate knowledge of the dynamic behavior of 
the spacecraft to ensure performance robustness to a 
variety of disturbances and uncertainties. 

In the past few years, several design approaches 
were proposed for vibration control during and after 
slewing maneuvers. NASA Langley Research Center 

initiated the Spacecraft Control Laboratory Experiment 
(SCOLE) program [V] to promote direct comparison and 
realistic test of various control design techniques 
against a common laboratory article. The article was 
intended to resemble a large space antenna attached to 
the space shuttle orbiter by a long flexible mast. 

The primary control objective of SCOLE is to direct 
the RF Line-of-Sight (LOS) of the antenna-like 
configuration towards a fixed target under the 

conditions of minimum time and limited control 
authority. 

The dynamical equations of slewing maneuvers of 
this large flexible spacecraft are developed by 
formulating Lagrange's equations using an inertial co- 
ordinate system and a body-fixed co-ordinate system at 
the point of attachment of the flexible beam to the 
shuttle. The generic model used for this analysis 
consists of a distributed parameter beam with two end 
masses. The three dimensional linear vibration analysis 
of this free-free beam model with end masses [2] is 
incorporated together with rigid-slewing maneuver 

dynamics [3] to yield the final set of highly nonlinear 
and coupled equations. 

The problem of directing the LOS of antenna-like 
configuration is being viewed as two phase control 
problem. In the first phase, the LOS of antenna-like 
structure is considered to be changed according to a 
prespecified target (slewing) together with 
minimization of vibration amplitudes to ensure 

stability. The second phase is to achieve the total 
vibration suppression at the end of slew maneuver with 
an augmented control law. 

The first phase of the control law is studied both 
in terms of pure rigid-body slewing and in terms of 
rigid-body slewing together with the suppression of the 
first two flexible modes. The dynamics of the motion 
during this phase are derived in terms of four Euler 
parameters, and using the method of nonlinear 
decoupling, this nonlinear slew maneuver control 
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problem Is reduced to calculating a pair of constants 
for each of the first three Euler parameters and the 
elastic modes in implementing output feedback. 

The second phase of the control is implemented at 
the end of the maneuver to achieve total vibration 
suppression, and this is performed by incorporating a 
state feedback control law derived from an infinite- 
time regulator problem formulation. 

Analytics: 

Kinetic Energy : 

The dynamics of slew maneuver are developed using 
two sets of co-ordinate systems and l.agrange's method. 
The body-fixed frame origin is located at the point of 
attachment of the flexible beam or mast (Fig. 1) . The 
second co-ordinate system is an inertial co-ordinate 
system. The transformation from the inertial fr'ame to 
the body-fixed frame is given by the matrix.C as 


cose 2 cos03 -cosegsines sine 2 


C= 


(sin 0 isine 2 cos 0 3 - {-sin 0 1 sin 0 2 sine 3 ' -sin0jCose 2 
4sin0^cos0,) +cos0^cos0j ) 

(-cos0isine 2 cos0 3 (coseisine 2 sine 3 cos0\cose2 


(1) 


Jj. sino 3 sinei) +cose 3 sin0i} J 

where if i,j,k represent the dextral set of orthogonal 
unit vectors fixed in the body-fixed fr^rne, then Q l is 
the rotation^jof f, 0 2 is the rotation of j and e 3 is the 
rotation of k. 

The angular velocity of the orbiter can be 
transformed from the inertial frame to the body-fixed 
frame for the body-three angles as 


W = M T e (2) 

The total kinetic energy expression of the system 
can be given as [4] 


T = T 0 + + T 2 


(3) 


where T 0 is the kinetic energy of the shuttle and is 
given as 

To= 1/2 mi V T V + 1/2 “> T i lU (4) 

The kinetic energy of the flexible beam is Tj and it 
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can be shown to be equal to 


T j * l/2mV T V n + 1/2 u T Ju- mV Jew* 1/2 d' T <Tdm 


+ V 0 T/d dm ♦uTftiS dm ♦ 1/2 [u x 'u *u ^ ]dl 

J "" rr. .1 


where 5 * c x and J is the mass moment of inertia 
matrix for the beam. The foregoing kinetic energy 
expression can be further simplified by using the 
three-dimensional modal analysis as 


Tj * l/2mV. 0 T V_+ 1/2 tu - mV_o^ ££ + m^ 9 1 - 2+ ^ T £ 


t* ... _ n *2 n •?- 
+ 0.1 + 1/4 p [ Z P q . + Z pg,qf ] 

i=l 5i i i=i 1 


where 


u x ■ 1 E 1 *xll s >1l< t > 


“* = i ! 1 +*i (s)q i (t) 


P 3 f <b . (s ) ds 


ii 6 


Pzi = / +vi( s ) ds 


P 3i = / S *xi (S) dS 


, 1 . 

p =J Sd> .(s) ds 

41 0 ; 


P 5i ’ / (s *xl » 2ds 
L 

P 61 = /( S * yj ) 2ds 

0 J 


I 


and 


a (t) 

lit) 


ijjl P 11 q 1 

i?l p 2i q i 
0 

“n . : 
i§l P 4i q i 

1$1 p 3i q i 


( 9 ) 


( 10 ) 


The kinetic energy T_ , of the tip mass (the reflector) 
is 2 

Tj = 1/2 m,vj Iq - BjvJltLJu ♦ nijvj d(L) 


- 1/2 m-s/H? (L)a(L)o * m, «Ta[L)d(L) 
t 1/2 m d_ T (L)i(L) + 1/2 n r I 2 n 


( 11 ) 


where 


R = <0 + 


“y '(L) 


(12) 


Equation (11) can be simplified as 
T = 1/2 m 2 V_ 0 T V_ 0 - m 2 vJt(L) w + m 2 vj dJL) 


l/2nfeJa(L)a_(L)(o + m 2 JajL)d(L) 


+ 1/2 
+ 1/2 P 1 I 2 P 


n n 


Ml jI}xi (L) *xj (L)q i q j t i Ii^ ( i L) *yJ (L)q i‘ , P 

P T !■> P + 1/2 d) T I jo (131 


where 
n 


n i 


( 13 ) 


( 14 ) 


iL=[.^ x1 ( L )q ,(t) .?j y (iL)q i(t) .^^L)q i (t)] 

Substituting T 0 » T i and \ from the foregoing equations 
into equation (3), the total kinetic energy expression 
can be written as 


T = 1/2 m 0 V_ T V_ + to 1 " HV + 1/2 io T I 0 w + V_ T A 1 g m 


+ “ T A iL + 1/2 i. A 3 q. 


(15) 


where 
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m 0 ■ m , ♦ P L + m 2 (16) 

H * ( p L ♦ n\j )jr + m 2 a( L) + P Lc ( 17) 

I 0 * I x * J ♦ I 2 (18) 

and also 


A 1 q » m 2 £(L) (19) 

A 2 1 ’ti + L + " 2r d J L ) ♦ n»2a(L)d(L) (20) • 


A 3 s 


\ u 
PL+m 2 + P 5i + P 6 i 


0 \ 

The matrix i T (L) is given as 


+ t T (L)I 2 *(L) 


(21) 



♦lx (L) ° ° 

0 *!y (L) 0 

0 0 * U < L > 


*U (L) 0 ° 

o *' 1y a) o 


( 22 ) 


L° 0 '♦#>] 

Equations of motion: 

Lagrange's equations for the case of independent 
generalized co-ordinates q k are 


1 ^ T - 9 T = Q k - 3 V (k = l, 2,....) 

dt FQ k 3 q k 8 Qk ( 23 ) 

where T(q,qJ is the kinetic energy, Vfql is the 
potential" energy and Qk are the generalized forces 
arising from nonconservative sources. 

However, the kinetic energy developed in equation 
(15) is given in terms of nonholonomic velocities V, u 
and generalized velocities C[. Using the chain-rule - 3 * 1/9 V 
be expressed in terms of generalized velocities. Also, 
if F(t) represents the total applied forces where 

F(t) =Io(t) +I 2 (t) (24) 
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then the generalized forces are given as CF(t). Thus, 
using the Lagrange's equations the translational 
equations can be obtained as 

d / 3T C T C ( U ) *f(t) (25) 

dt v ay ; 

which can be simplified as 

n^V_- H i + A^'fi + fCt) (26) 

where the nonlinear term Njis given as 

N i « -C T C ( m oJ ” H w+ A i q ) (27) 

= “ (m oX." H « + Aijq ) 

Similarly, using equation (2) and the chain rule in the 
Lagrange’s equations, the rotational equations are 
obtained as 

HV + I 0 y+ A 2 q = G(t) + N 2 (28) 

where Gft) is the net moment about the mass center of 
the orEiter and is given as 

G = G + ( r + a ) x F (29) 

_ _ o - — -2 ' ' 

and the nonlinear term is given in terms of 

transformations M and C, and y ,1 and £. .The vibration 
equations of the beam can be obtained by again using 
Lagrange's equations and the potential energy function 

U = 1/2 q T l< q (30) 


where the stiffness matrix K is given as 
6 i , 4 


X = 


UX.1U 
L 3 


The vibration equations are 
A + A^y. + A 3 q = -Kq 


(31) 

(32) 


Slew Maneuver Dynamics: Considering the translational 

velocity and acceleration to be negligible during the 
slew maneuver, the rotational equations and vibration 
equations are 

Io £L + A 2 £ = _G(t) + Ji 2 ( y) 

T. 

As “ + A 3 q = -Kq 
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(33) 

(34) 


( 35 ) 


Equation (33) can be rewritten as 

£ = Io“ 1 C £ + JL *( “) “ A 2 q 3 

The first three Euler parameters are defined as 


ci 
e 2 

L e 3 J 


= x. sin i>/2 


ei 4 = cos i>/2 

d£ = 1/2 ( G 4 10 + £ x <d) 
dt 

de^ 55 "1/2 s 

dt , , dc . de . 

d£ dt 

i = St = iL ( £» “) 

Defining the output vector as 


L = L 

y_ = £ = h_( £ , £ ) 

• • 3 h • 3h • 

y = _ — £ + £ 

~ 3e_ 3£ 

y = M £ . £ ) + p ( £ . “ .q_ )<L 

Choose G as 

G = P -1 {- £+ £ } 

It can be shown that ^always exists. 
Then 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 


1*1 =t (47) 

form a system of uncoupled equations and they can be 
expressed in terms of output feedback as 

y_= [k p ]y_ + [k v ] l (48) 

The elements of these gain matrices can be chosen for 
desired system response. 

Results 0 

In figures 2 and 3 , the slewing maneuver of 20 is 
shown in roll plane. The maneuver shown in figure 2 is 
for a pure rigid-body slew maneuver whereas figure 3 is 
the same maneuver with suppression of two first 
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flexible modes. The corresponding case in the pitch 
plane is shown in figures 4 and 5. A rigid-body slew 
maneuver is shown about an arbitrary axis in the 
remaining figures. The effect of first four flexible 
modes are incorporated in this slew maneuver and 
figures 6, 7 and 8 in this case represent the moment 
components to perform this maneuver. 

Conclusions 

The equations of motion of SCOLE model are highly 
nonlinear and coupled and this results into the 
excitation of higher flexible modes when the lower 
modes are being controlled during the slewing maneuver. 
Although for this study the beam vibrations at the end 
of slew maneuver are controlled using Linear Infinite- 
time Regulator Problem formulation, analysis of higher 
uncontrolled modes indicates serious control spillover 
problem due to coupling among the modes. 
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Fig. 4 Slew Angle vs. Time in Pitch Plane 

(Rigid -Body Model) 
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Slew Angle in-D 



: ig. 5 Slew Angle vs. Time in Pitch Plane 
(Model with Two Flexible Modes) 
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Moment Component G 3 vs. Time 
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